LONG-TIME ASYMPTOTIC FOR THE DERIVATIVE 
NONLINEAR SCHRODINGER EQUATION WITH 
DECAYING INITIAL VALUE 



JIAN XU AND ENGUI FAN* 

Abstract. We present a new Riemann-Hilbert problem formal- 
ism for the initial value problem for the derivative nonlinear Schrodinger 
(DNLS) equation: 

iq t (x,t) + q xx (x,t) + i{\q\ 2 q) x = 

on the line. We show that the solution of this initial value problem 
can be obtained from the solution of some associated Riemann- 
Hilbert problem. This new Riemann-Hilbert problem for the DNLS 
equation will lead us to use nonlinear steepest-descent/stationary 
phase method or Deift-Zhou method to derive the long-time as- 
ymptotic for the DNLS equation on the line. 



1. Introduction 

The main purpose of this paper is to develop an inverse scattering ap- 
proach, based on an appropriate Riemann-Hilbert problem formulation, 
for the initial value problem for the derivative nonlinear Schrodinger 
(DNLS) equation [9J on the line, whose form is: 

iq t {x,t) + q xx {x,t) - i(rq 2 ) x = 0, (1.1a) 

q(x,0)=q (x). (1.1b) 

where r = ±q, q denotes complex conjugate of q, the subscripts denote 
differentiation with respect to the corresponding variables. And in this 
paper we use r = —q 
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The DNLS equation has several applications in plasma physics. In 
plasma physcis, it is a model for Alfven waves propagating parallel 
to the ambient magnetic field, q being the transverse magnetic field 
perturbation and x and t being space and time coordinates, respectively 

DO- 

Our goal is to develop the inverse scattering approach to the DNLS 
equation, in view of its further application for studying the long-time 
asymptotics. The starting point of the approach is the Lax pair repre- 
sentation: the DNLS equation is indeed the compatibility condition of 
two linear equations [9]: 

v i x + ik 2 v\ = qkv2, v 2x — ik 2 V2 = rkvi, (1.2a) 

iv u = Avi + Bv 2 , iv 2t = Cvi - Av 2 , (1.2b) 
where £ G C is the spectral parameter, and 

A = 2k 4 + k 2 rq, B = 2ik 3 q — kq x + ikrq 2 , 

C = 2ik 3 r — kr x + ikr 2 q. 

In the present paper, we propose a scattering inverse scattering for- 
malism, in which the Lax pair is used in the form of a system of first 
order matrix-valued linear equations. Then dedicated solutions of this 
system are defined and used to construct a Riemann-Hilbert (RH) prob- 
lem in the complex plane. The main advantage of the representation of 
a solution of the DNLS equation in terms of the solution of a Riemann- 
Hilbert problem is that it allows applying the nonlinear steepest descent 
method by Deift and Zhou [6] in order to obtain rigorous results on the 
long-time asymptotic behavior of the solution. 

An alternative inverse scattering method based on a Riemann-Hilbert 
problem formulation for the DNLS equation can be founded in [T9J for 
the Cauchy problem, and in [1] for the initial-boundary value problem 
on the half-line, in [23] for the initial-boundary value problem on the 
interval. 

In Section 2, we define appropriate eigenfunctions and spectral func- 
tions, which are used in Section 3 in the reformulation of the scattering 
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problem as a Riemann-Hilbert problem of analytic conjugation in the 
complex plane of the spectral parameter. And in Section 4 we obtain 
the long-time asymptotic behavior of the solution of DNLS equation 
by the method of the Deift-Zhou/nonlinear steepest descent based on 
the new Riemann-Hilbert problem which is obtained in the subsection 

2. ElGENFUNCTIONS AND SPECTRAL FUNCTIONS 

2.1. Eigenfunctions. First introducing 

*-(:)■ 0-0- 

we can rewrite the Lax pair fll.2p in a matrix form: 

i/j x + ik 2 a 3 i/j = kQijj, 

i/H + 2ik 4 a 3 i{> = {-tk 2 Q 2 a 3 + 2k 3 Q - ikQ x a 3 + kQ 3 )^j, 
Extending the column vector ip to a 2 x 2 matrix and letting 

^ = , l p e i(k 2 x+2kH)a 3 

we obtain the equivalent Lax pair 
V x + ik 2 [a 3l ty] = kQ^, 

% + 2ik 4 {a 3 , *] = (-ik 2 Q 2 a 3 + 2k 3 Q - ikQ x a 3 + kQ 3 )^ {2.2) 
which can be written in full derivative form 

d(e i{k2x+2k4t)f7:i ^(x, t, k)) = e l{k2x+2kH) ^U(x, t, jfc)*, (2.3) 

where 

U = U 1 dx+U 2 dt = kQdx+(-ik 2 Q 2 a 3 +2k 3 Q-ikQ x a 3 +kQ 3 )dt. (2.4) 

In order to formulate a Riemann-Hilbert problem for the solution of 
the inverse spectral problem,we seek solutions of the spectral problem 
which approach the 2x2 identity matrix as k — > oo. It turns out that 
solutions of Eq.f l2.3p do not exhibit this property [9], hence we have 
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to transform the solution \l> of Eq.f l2.3p into the desired asympototic 
behavior [3]. 

Consider a solution of Eg. ( 12. 3 p of the form 

\Pl \&9 ^/q 1 

9 = D + T + 'W + 'i + 0{ ¥ )t k ^°° 

where D, ty 2 , ^3 are independent of k. Substituting the above ex- 
pansion into the the first equation of (I2.2p .and comparing the same 
order of fc's frequency, it follows from the 0(k 2 ) terms that D is a di- 
agonal matrix. Furthermore, one finds the following equations for the 
0(k) and the diagonal part of the 0(1) terms 

0{k):i[a 3 ,* 1 ] = QD, i.e. ^ = l -QDa 3 , 

with being the off-diagonal part of \&i,and 

0(1): D x = Q¥°\ 

i.e. 

D x = % -Q 2 o 3 D. (2.5) 

On the other hand, substituting the above expansion into the second 
equation of (12.21) . one obtains from that 

0(k 3 ):2t[a 3 ^ 1 } = 2QD, i.e. #< o) = % -QDo 3 - (2.6) 

and 

0{k) : 2i[a 3 , * 3 ] = -iQ 2 <j 3 ^ ] + 2QVif - iQ x a 3 D + Q 3 D, (2.7) 

i.e. 

- iQ 2 a^f + 2Q¥ 3 o) = —Q 3 ^ + ~QQ X D + % -Q A o 3 D, (2.8) 

where ^ denotes the diagonal part of and for the diagonal part 
of the 0(1) terms 

0(1) : D t = -iQ 2 ^ + 2Q¥° ] - iQ x a 3 ^ + Q 3 ^ , 

again, using ( 12. 6 p and ( 12. 8 p . we have 

3? 1 
D t = (-Q i a 3 + -[Q,Q X ])D, 
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which can be written in terms of q and r as 

3i 1 

A = ( jr 2 g 2 + -(r x q - rq x ))a 3 D. (2.9) 

We note that Eg. ( 11.1 all admits the conservation law 

(^)t = (jr 2 q 2 + ^(r x q - rq x )) x . 

Because we just consider the Cauchy problem for the DNLS eguation 
fll. lap , the two Egs.( l2.5p and (I2.9P for D are consistent and are both 
satisfied if we define 

D(x,t) =e lJ (-^ Aa \ (2.10) 

where A is the closed real-valued one-form 

1 3 i 

A(M) = -zqrdx + (~r 2 q 2 - ~(r x q - rq x ))dt. (2.11) 

Noting that the integral in (12. lOf) is independent of the path of integra- 
tion and the A is independent of k, then we introduce a new function 
H by 

ty(x,t,k) = e if (+£*) Aa3 fi(x,t, k)D(x,t), (2.12) 

Thus,we have 

u = l+0( T ), k^oo, (2.13) 
k 

and the Lax pair of Eq. fl2.3j) becomes 

d{e i{k2x+2kH)a3 fi{x, t, k)) = W{x, t, k), (2.14) 

where 

W(x, t, k) = e ^ 2 x+2kh)a A y^ x ^ tj ^ 

with 

V = Vxdx + V 2 dt = e~ iJ (+«.*3 A&3 (U - iAa 3 ). 
Taking into account the definition of U and A, we find that 

(x,t) 



--rq kqe J (+°°>t) 



kre J (+°°.*) ~rq 
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—ik 2 rq — ^r 2 q 2 — \[r x q — rq x ) (2k 3 q + ikq x + kq 2 r)e 2l ^(+ 



(x,t) 
-oo,t) 



(2k 3 r + ikr x + kr 2 q)e 2l ^+ ao - t '> A ik 2 rq + ^r 2 q 2 + \{r x q — rq a 

(2.16) 



Then Eq. fl2.14p for fi can be written as 

fj. x + ik 2 [a 3 ,/j] = Vxn, 
fit + 2ik 4 [a 3 ,fi] = V 2 fi. 



(2.17) 



Throughout this section we assume that q(x, t) is sufficiently smooth, we 
define two solutions of Eq. ( I2.14p by 



Mx,t) 

fj,j(x, t,k) —1+ e 



- ,,kx+2kt)&3 W(y,T,k), J = 1,2, (2.18) 
where (xi,ti) = (+00,^,(^2,^2) = (— 00, t),see Figure 1. 




O 



Figure 1. Paths integrals of the fii and fi 2 - 

The analytic properties of the 2x2 matrices fij(x,t;k), j = 1,2, 
that follow from (I2.18p . are collected in the following proposition. We 



denote by /ij\x,t, k) and ^\x,t,k) the first and second columns of 
fij(x,t; k), respectively. 

Proposition 2.1. The matrices fii(x,t; k) and ^{Xjt; k) have the fol- 
loeing properties: 

(i) det /ii(x, t, k) = ^{x, t; k) = 1. 

(ii) jj!~l\x,t,k) is analytic in \m.k 2 < and 



,( 2 ) 



1 I + 0{\), as k ->■ 00, ImA; 2 < 0. 
k 
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(iii) ^ {x,t,k) is analytic in Irak 2 > and 

fji?\x,t,k) = ( ° ) +0(t), as k ->■ oo, ImA; 2 > 0. 
V 1 / fc 

(iv) /j,2 (x,t,k) is analytic in Irak 2 > 0, and 
H { 2 1] {x,t,k) = ( J j +0(i),as fc -»■ oo, ImA; 2 > 0. 

(v) fi^\x,t,k) is analytic in Irak 2 < 0, and 
H^\x,t, k) = [ J J +0(i),as fc ->■ oo, ImA; 2 < 0. 



(vi) Moreover, 



as k — y oo a/ong curves transversal to the real and image axis, 
where 



[a 3 ,ji(x,t)] 

and 



q(x,t) 



-q(x,t) 
q(x,t) = q(x,t)e' 2lJ +^ A (2.19) 



Since the eigenfunctions fi\(x,t, k) and /i 2 (x,t, k) satisfy both equa- 
tions of the Lax pair, we have 

fi 2 {x, t, k) = fii(x, t, k)S(k), k 2 eR, (2.20) 

where S(k) is independent of (x, t) and is defined in (12. 23j) 

Proposition 2.2. (Symmetries) For j = 1,2, the function fi(x,t,k) = 
Hj(x,t,k) satisfies the symmetry relations: 



fj. u (x,t, k) = fi 22 (x,t,k) , ^ 22 ^ 
H2i(x,t, k) = ni 2 (x,t, k), 
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as well as 



fJ>2l{X, t, - 



k) 
k) 



k) 



k) 



fM n (x,t,k), 
-fi 12 (x,t, k) 
-fJai(x,t, k) 
Vmfat, k). 



(2.22) 



Thus, we have 



S(k) 



a(k) -b{k) 



b{k) a(k) 



) 



(2.23) 



3. The basic Riemann-Hilbert problem 



3.1. The original Riemann-Hilbert problem. The scattering re- 
lation ( 12 .2 Op involving the eigenfunctions "$?(x, t, k) = fii(x,t,k) and 
<&(x,t,k) = fi2{x,t,k) can be rewritten in the form of conjugation of 
boundary values of a piecewise analytic matrix-value function on a 
contour in the complex k— plane, namely: 



where M±(x, t, k) denote the boundary vales of M(x, t, k) according to 
a chosen orientation of S, and S = K U zR 

Im(k) 



M. 



+ (x, t, k) = M_(x, t, k)J(x, t, k), 



k 2 E R, 



(3.1) 



Re(k) 



Figure 2. The jump contour for M. 
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Indeed,let us write (12.201) in the vector form: 



a(k) 
& 2 \x,t,k) 
a(k) 



^ {1 \x,t, k) + r(k)¥ 2 \x,t,k), 



(3.2) 



-r(k)^ {1) {x,t, k) + ¥ 2) {x,t,k), 



where 



and define the matrix M(x, t, k) as follows: 

( ( * (1 Hy*) e «g(fc) mV)(x,t,k)e- it6{ V ), k e {k e C\lmk 2 > o}, 

M ^ k) = | ( *V)( x ,t,k)e«°W ^Mhle-^ ), ke{ke C\lmk 2 < 0}, 



(3.4) 



where 



9(k) := 2k 4 + -k 2 , (3.5) 

t 

Then the boundary values M+(x, t, k) and M_(x, i, k) relative to E are 
related by (13. ip . where 



, l-r(k)r(k) - r (k)e- 2iW ^ . 
■/(*,*,*)= I r(Jfe) ^ ( V l \ I. (3-6) 

The jump relation (13. ip considered together with the properties of 
the eigenfunctions listed in Proposition 12 . 1 1 suggests a way of represent- 
ing the solution to the Cauchy problem ( II. ip in terms of the solution 
of the Riemann-Hilbert problem, which is specified by the initial con- 
ditions (ll.lbp via the associated spectral function r{k). 

The function q(x,t) can be expressed in terms of the solution of the 
basic Riemann-Hilbert problem as follows: 

q(x, t) = 2i lim (kM(x, t, k)) n . (3.7) 

where M is the solution of the following Riemann-Hilbert problem: 
The original Riemann-Hilbert problem: Given r(k),k 2 G R, and 
S = 1U z'R, find a 2 x 2 matrix- value function M(x, t, k) such that 
(i) M(x,t, k) is analytic in k e C\S. 
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(ii) The boundary value M± (x, t, k) at S satisfy the jump condition 

M + (x,t,k) = M-(x,t,k)J(x,t,k), fceS, 

where the jump matrix J(x, t, k) is defined in terms of r(k) by 
dSU). 

(iii) Behavior at oo 

M(x,t,k) = I + 0(~), as k -» oo. 

And from the definition of the function q(x, t) in (I2.19P we find 

\q\ = \q\ (3.8) 

this means that the solution the the initial value problem (11.11) can be 
expressed as follows: 

q(x, t) = q(x, t)e~ 2i I«(**)l a «& (3.9) 

3.2. The new Riemann-Hilbert problem. The jump condition (13. 6p 
is obtained in [19]. In that paper, the authors used this condition (13.61) 
to analysis the long-time asymptotic behavior. But if we try to analysis 
the long-time asymptotic behavior of the DNLS equation (ll.lal) with 
step-like initial value problem, this type of Riemann-Hilbert problem 
has a contradiction in the plane wave region. So we try to derive a new 
Riemann-Hilbert problem, which is similar to the type of nonlinear 
Schrodinger equation, to overcome this contradiction. In this paper, 
we just analysis the long-time asymptotic behavior of the DNLS equa- 
tion with decaying initial value problem. The step-like initial value 
problem will be obtained in another paper. 
We define 

N(x, t, k) = k~^M(x, t, k), (3.10) 
then the jump condition for iV is 

,-i(k?x+2kH)a> ( 1 - r(k)r(k) 



N+(x,t,k) = N-(x,t, k)e 



kr(k) 
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introducing A = k 2 and control the branch of k as Sign \mk =Sign ImA, 
and define the modified scattering data p(A) = [9]. And defining 

X = DX . where D = I 1 ° I and b = £ 



6 1 



2; : 



A 

Figure 3. The jump contour for N. 
Then the jump condition for X is 

X+(x,t,X) = X_(x,t, \)e~ l(Xx+2x2t)&3 J N {x,t, \). (3.12) 

where 

/ l-A|p(A)| 2 -p(\) 
J <^ ( ' A >=( XfAX) ! 

the matrix Jjy admits the following triangular factorizations: 



1 




-p 



and the solution of the DNLS equation (jl.lj) is 



1-A|p|* 

1 



A— >oc 



(3.13) 



q(x, t) = 2i lim (AA^(x, t, \)) 12 , (3.14a) 



q(x, t) = q(x, t)e- 2i l^*)l 2 *. (3. 14b) 

4. Long-time analysis 

In order to analysis the long-time behavior of the solution of the 
DNLS equation, first we should split the jump matrix into an appro- 
priate upper /lower triangular form, then this can help us localize the 
problem to the neighborhood of the stationary point. An appropriate 
rescaling then reduces the problem to a Riemann-Hilbert problem with 
constant coefficients, which can be solved explicitly in terms of classical 
functions. 
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4.1. The first transformation. In this subsection we extend the 
Riemann-Hilbert problem ( 13 . 1 2 j) to an augmented contour of the type 
Si given in Figure 4, which is constructed to reflect the signature of 
R,eit8. For technical reasons we will assume that the contour Ex is 
composed of straight lines, as shown in Figure 4 below with angle j, 
although any contour of the same general shape as Ei would do. 



Figure 4. The contour of Ei. 

As in the [6], we first consider the stationary point of the function 
0(A) = 2A 2 + f A, that is, letting 

dX 

we get the stationary point A = = rr. 

And we also get the signature table of Re#(A) that is as follows in 
Fig.5. 

Then we introduce a scalar function. Let 5(A) be the solution of the 
scalar factorization problem 

<5 + (A) = MA)(l-A|p(A)| 2 ), A<A , 
5 + (A) = MA), A>A , (4.1) 

5(A) -> 1, A -)■ oo. 

Direct calculation shows that (14.11) is solved by the formula 

, (A) * r° '°g(i-/if>iv (4. 2) 

2m J_ 00 A' - A 

And we can find that 5(A) and 5(A) _1 are uniformly bounded in A 
and for |Aq| < M. 
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Re{i6) > 



Re(i9) < 



Re{i6) < 



Re(X) 



Re(i6) > 



Figure 5. The signature table of Re(z0(A)). 



We conjugate the Riemann-Hilbert problem (I3.12p by 

r>W ' S(X) 



5-\X) 

leads to the factorization problem N^(x,t, A) = N(x,t, A)5~°" 3 (A), 



(4.3) 



N ( +\x,t,X) = (x, t, A) J N (i) (x, t, A) , AgI, 



(4.4) 



where 



1 



1-A|p|* 



J N (i)(x,t, A) = 5™J NW (x,t, A)5; CT3 (A) = < 



-r2 



1 -pS 



1 



1 

\ P 5~ 2 1 
(4.5) 

Having made the above definitions, we now describe the strategy. 
Suppose that the coefficients 



A < A 
A > A 



[P], [Ap] (4.7) 
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can be replaced by some rational functions 

[ 1-A|p| 2j ' [ 1-A|p| 2j ' 
respectively. Then if the poles of these functions are appropriately 
placed, the Riemann-Hilbert problem on R can be deformed to the 
contour Si. 

Remark 4.1. In this paper, we assume that the function p(A) has no 
zero. 

To verify that the coefficients ( 14. 6 p can be replaced by the rational 
functions ( 14. 7p with well-controlled errors, we proceed as follows. 
1: For A < A =2 



°' 1-A|p| 2 

Set 



1 - X\p 

By Taylor's formula, we have 



( A -i)-+5/(A) = yUo+ ^ 1 (A-A )+- • .+^ m (A-A ) m +i- / A ((.-^)-+5 / (-)) (m+1) (7)(A-7) m rf7 

m! ./),„ 



Ao 

(4.9) 



and define 



fc(A) = /(A)-i2(A), (4.11) 
As before, the proof of the following result is straightforward: 



Lemma 4.2. 

*'/(A) _ d?R{\) 
d\i |a ° ~ dAi |a ° 
t4/so, yUj = /ij(A ) decays rapidly as A — >■ oo 



< j < m. (4.12) 



Proof. Formula (14.121) is immediate. The decay as A — > oo 
follows from the formulae 

^ i! tZti* |Ao ' 

□ 
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In what follows we fix m 6 Z + and, for convenience, we 
assume that m is of the form 

k = 4q + l, qeZ + , (4.13) 

Write 

f(\) = h(\)+R(\), A<A , (4.14) 
Then by the above lemma, 
d j h(X) . 

-^Uo = 0, 0<j<m, (4.15) 

We now use this property to split h further as 

h(X) = h 1 (X) + h 2 (X), (4.16) 

where hi(X) is small and /12(A) has an analytic continuation to 
A + iO. Thus 

f = h 1 + (h 2 + R) (4.17) 

is the desired splitting of /. 

Set 

Consider the Fourier transform with respect to 9. As A — > 9(X) 
is one-to-one in A < A ,we define 

= 0, 9 < -2X 2 

Thus, as 9 > — 2Ag, from formulae (O) . fl4TTTj) and fl4TT8|) it 
follows that 

| (A)= (A-i)^ g(A - Ao) - (4 ' 20) 

where 

,Ao) = ^r [\(--i) m+5 f(-)Y m+1 \X +u(X-X ))(l-urdu. (4.21) 



from which we see that 
d 3 g{X, A ) 



< C, A < A . (4.22) 
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I. 



00 

oo 



Then, we obtain 

2 



dd 



Xo 



r 

< c i- 



'(|(A)) 4(A-A )dA 



A 



(A-Aq) 



771+1 — 5 — 2,7' 



(A-i) m + 3 -9 



(A - A )rfA < C 



f or o < 7 < ^±i^ = 3i±2. 



By Plancherel, 



l + s 2 ) j \(h/f3){s)\ 2 ds < C < oo, 0<j< 



3g + 2 



A„ 



ise(A) (V/3)(A)^(A), 



where 
And by Fourier, 

/■oo 

(h/P)(\)= / e isew {h/P)(s)ds. 



In the above formulae we use the convenient notation <is 
and d(6(X)) = ^ 



(4.23) 
(4.24) 

(4.25) 

(4.26) 

ds 
/2n 



Remark 4.3. The constants in (i4-23\) and ( fj . 24\ ) should prop- 
erly be denoted by C\ and c 2 . Here, and in what follows, we use 
c and sometimes C to denote a generic constant. This abuse 
of notation should not give rise to any confusion. 

We split 

h(\) = /3(A) J t °° e is6 ^{h/P)(s)ds + /3(A) e is0 ^{hj%{s)ds 



= h 1 (X) + h 2 (X). 
For A < Aq we find that 



- 2it6 ^h x {X)\ < 



C 



for any p < 



3q + 2 



(4.27) 



(4.28) 



\x-i\ 2 tp-r " ' 2 

On the other hand, /12(A) has an analytic continuation to 
the upper half-plane, where Rez6>(A) > 0, and for A on the line 
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■ 37T 

A + \ ue l ~ , u > 0, 

\e- 2MW h 2 (X)\< C J x e _. lg+2 . (4.29) 
However, from expression of 0(X), that is 

0(A) = 2(A - A ) 2 - 2\l (4-30) 

we have 

Rei9(X) = 2X 2 u 2 , (4.31) 

and hence 



\ e -2ite(\) h mi < c\l\{(t\l)huYe-^ 2 ] 

1 - |A-i|«+=»(*Ag)* ( 4>32 ) 



< s < _ 

~ \\-i\i+ 2 ti ~ \\-i\ 2 ti ' 

On the line Aq + Xoue 1 ^ , u > e, e > we have 



\ e -x*WR(X)\ < Ce- 4tx o u2 < Ce^ tx o £2 (4.33) 



2: For case: A < A , fz^p 
Set 



AW = j^p. (4.34) 



Again, by Taylor's formula, we have 

(A+0 m+5 /(A) = ^o+/ii(A-A )+- • •+ / u m (A-Ao) m +^ T / A ((-+ ? ) m+5 /(-)) {m+1) (7)(A-7) m rf 7 , 

(4.35) 

and define 

h(\) = f(\) - R(\), (4.37) 
As before, the proof of the following result is straightforward: 



Lemma 4.4. 

^#(A) 

Also, ^ = fJ>i(Xo) decays rapidly as A — > oo. 



A , 0<j<m. (4.38) 
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Clearly ^B| Ao = 0,0 < j < m. 



Set 



/3(A) 



(A - A )< 



(A + i)?+2- 

From formulae fl4~35l) . (07j) and (14391) it follows that 

(A-A r +i -« 

where 



(4.39) 



(4.40) 



^(A,A ) = ^ / ((■+^r +5 /(-)) (m+1) (Ao+n(A-Ao))(l- M ) m ^. (4.41) 
ml 



from which we see that 
d j g(X,X Q ) 



<C, A < A . 



(4.42) 



Then, we obtain 

2 



Ao 

oo 



(A-Ap) 



m+l — q-2j 



(A+i) m + 3 -"J 

f or o < j < = 3^2. 



(A - A )rfA < C 
(4.43) 



By Plancherel, 



(l + s 2 ) j \{h/P){s)\ 2 ds<C <oo, 0<j< 



3g + 2 



where 

And by Fourier, 

(V/?)(A) 
Again we split 



Ao 



> e(A) (V/3)(A)^(A), 



- is9 ^(h/p){s)ds. 



(4.44) 



(4.45) 



(4.46) 



h(\) = /3(A) /°° e- iae W(h/ '/3)(s)ds + /3(A) f^e-^ih/ '/3)(s)di 



h 1 (X) + h 2 (X). 



(4.47) 
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For A < Ao we find that 

\e^ hl (X)\ < - C 1 ; for any p < **±1. (4.48) 

\X + l\ 2 t p 2 z 

On the other hand, h 2 (A) has an analytic continuation to 
the lower half-plane, where Rei#(A) < 0, and for A on the line 

A + X ue~ l ~ , u > 0, 

\e^h 2 (X)\ < C j" i{q+2 ■ (4.49) 
However, from expression of 9(X), that is 

0(A) = 2(A - Ao) 2 - 2A 2 , (4.50) 

we have 

Rei6(\) = -2X 2 u 2 , (4.51) 

and hence 

\ e 2ite(\) h / a ni < c\l[{{t\l)huYe- 2tx h 2 ] 

\\+i\«+ 2 {t\l)i (4.52) 

< C -^r < 



\x+i\i+m — \\+i\ 2 ti 



: 3tt 



On the line Ao + Ao«e '4 ) «>e,£>0we have 

\e- 2U9 ^R(X) \ < Ce- 4tx o u2 < Ce- Atx ^ 2 (4.53) 

In fact this case is just the conjugate of the above case. And 
the two cases in the following is fimilar with these two cases, 
but we write them down here for the reader's convenience. 



3: For case: A > A , p(A) 
Set 



/ 3 (A) = p(A), A>A , (4.54) 
Again, by Taylor's formula, we have 

(A _, r +5 /(A) = /Wi(A _ Ao) + . . . +fXm{X _ Xor+ 1 /■ A (( ._ i) m + 5 /( . ))(m+ l) (7)(A _ 7r ^ 

m - JXo 

(4.55) 

and define 
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h(\) = f(\)-R(\), (4.57) 
As before, the proof of the following result is straightforward: 



Lemma 4.5. 



d j f{\) . d j R(X) 



I Aq 3 



< j < m. 



(4.58) 



Also, Hi = /ij(A ) decays rapidly as A — >■ oo. 



Clearly 
Set 



A = 0, < j < m. 



(A-i)«+ 2 ' 

Prom formulae ( 1435]) . (T4~57]) and (T4~59]) it follows that 



(4.59) 



(4.60) 



where 



g(K A c 



((-O m+5 /(-)) (m+1) (Ao+«(A-A ))(l-n) m rf M . (4.61) 



from which we see that 
^^(A, A ) 



oo 
oo 



4lW*(A(0))) d6= J* 



<C, A > A . 



(4.62) 



Then, we obtain 

2 



^AV7A(A)) 4(A-A )JA 



-4(A-Ao)dA/ V/8 1 
(A-Ao) m+1 -g- 2 J 



(A-i) m + 3 "'3 



for < j < = ^±1 

By Plancherel, 

o 

(l + s 2 y\(h//3)(s)\ 2 ds <C< oo, 0<j< 

X) 

where 



(A - A )rfA < C 
(4.63) 

(4.64) 



(tyfl( a ) 



is0W {h/p){\)d6{\), 



(4.65) 



A 
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And by Fourier, 

/oo 
e isd ^\h/l3){s)ds. (4.66) 
■oo 

Again we split 

h(X) = /3(A) J t °° e is9 ^{hjf3)(s)ds + /3(A) e ise ^{h/~f3)(s)ds 
= h l (X) + h 2 (X). 

(4.67) 

For A > A we find that 

|e- M * (A) /n(A)| < — — — j-, for any p < ^L±l, ( 4 . 68 ) 

|A — i\H p 2 2 

On the other hand, /12(A) has an analytic continuation to 
the lower half-plane, where Kei6(X) > 0, and for A on the line 

A + X ue~ l %, u > 0, 

\e-^h 2 (X)\ < C J x _ i]q+2 • (4.69) 
However, from expression of 0(A), that is 

0(A) = 2(A - Ac) 2 - 2A 2 , (4.70) 

we have 

Rei6(X) = 2X 2 u 2 , (4.71) 

and hence 

\ e -2ite(\) h (x) \ < c\l[{(t\l)huYe-^ 2 ] 

\*-i\«+ 2 (t\ 2 )i (4.72) 

< £ < c 

~ ~ \\~i\ 2 ti ' 

On the line A + X ue~ l ^,u > e, e > we have 

\e- 2ite ^R(X)\ < Ce- 4U o" 2 < Ce'^ (4.73) 

4: For case: A > A : Ap(A) 
Set 

/(A) = Ap(A). (4.74) 
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Again, by Taylor's formula, we have 



(A+2) fc+5 /(A) = /io+MA-Ao)+- ■ -+/i m (A-Ao) fe +- 



ml 



and define 



R(X) 



(A + *) m + 5 ' 
h(X) = f(X)-R(X), (4.77) 
As before, the proof of the following result is straightforward: 



((■+^ +5 /(-)) (m+1) (7)(A-7r^7, 
(4.75) 

(4.76) 



Lemma 4.6. 



<Pf{\) . d^R(X) 



IAq 



I Ao 3 



< j < m. 



(4.78) 



yMso, /ij = /ij(Ao) decays rapidly as Ao — >■ oo. 
Clearly ^B| Ao = 0,0 < j <m. 



Set 



/3(A) 



(A - A )'- 



(A + i)o+ 2 

From formulae fl4~T5l) . fl4~TTj) and ffl~79]) it follows that 
7I (A) - (A + z)™+3-^ (A ' Ao) ' 



(4.79) 



(4.80) 



where 



#(A, A c 



m! 



'((.+i) m+5 /(-)) (m+1) (Ao-H*(A-Ao))(l-«) m dti. (4.81) 



from which we see that 
d?g{\, A ) 



< C, A > A . 



(4.82) 



Then, we obtain 



(A-A ) m+1 - g - 2j ' 



(A+«) m+3 ~' J 



(A - A )rfA < C 



for < j < = ^±1. 



(4.83) 
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By Plancherel, 

r \l + s 2 ) j \{h/^)(s)\ 2 ds <C <oo, 0<j<?H±l. (4.84) 

J-oc 2 

where 

/>oo 

(/i//9)(s) = / e ise W(h/(3)(\)d9(\), (4.85) 
And by Fourier, 

/oo 
e -MW(h/p)(s)ds. (4.86) 
-oo 

Again we split 

/i(A) = /3(A) / t °° e- ise ^(h/p)(s)ds + /3(A) e" is 9 ^(hjf3)(s)ds 
= h 1 (X) + h 2 (X). 

(4.87) 

For A > Ao we find that 

|e 2 ^ A )/n(A)| < — — — y, for any p < - + 1 ~ 9 . (4.88) 

|A + i| 2 t p 2 2 

On the other hand, /12(A) has an analytic continuation to 
the upper half-plane, where R,ei6(X) < 0, and for A on the line 

A + Xoue 1 ^ , u > 0, 

\e^h 2 (X)\ < ,°" + e . |g+2 • (4.89) 
However, from expression of 0(A), that is 

0(A) = 2(A - Ao) 2 - 2A 2 , (4.90) 

we have 

Rei6(X) = -2X 2 u 2 , (4.91) 

and hence 



\ e 2ite(X) h < cAg[((tAg)? u )g e - 2 ^Q" ] 

< £ < £ 

— |A+i|«+ 2 t$ ~~ |A+i| 2 <2 ' 

On the line A + X ue l % , u > e, e > we have 

\e- 2itew R(X) \ < Ce- 4U o" 2 < Ce- 4U o £2 (4.93) 
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We can summarize our results as follows: let / be an arbitrary positive 
integer and let k — 4g + 1 be sufficiently large that the integers that 
are the last formula of the above formulas about h 2 are all greater than 
/. Let L denote the contour 

L : {A = A + X ue-^;u > 0} U {A = A + Xoue^; u > 0} (4.94) 

so that the contour Ei in Figure 4 is given by 

£i = LULUR. (4.95) 

Also set 

L £ = { A = A + Aowe'T , w > £:} U { A = A + X ue~^ ,u>e} (4.96) 
Proposition 4.7. Let 



/(A) 



/2(A) 



Xp 



1-X\p\- 



-, A < A c 



/ 3 (A)= p(A), A>A 
{ / 4 (A) = Ap(A), A > Ac 

Tnen / /tas a decomposition 

f(X) = h 1 (X) + h 2 (X) + R(X), A el 



(4.97) 



(4.98) 



where R(X) is piecewise rational and /12(A) aas an analytic continuation 
to L or L satisfying for case &4-l\) and 



-* tf WAi(A)| < 
- 2 ^W/i 2 (A)| < 



(1 + |A| 2 )t r 



and 



-2it6{\) 



R(X)\ < Ce 



or /or case OT71) ana 1 OT7 



3 2 ^W/n(A)| < 

3 2U0 ^h 2 (X)\ < 



(1 + |A| 2 )^' 



(l + |A| 2 )t r 



A g R, 
A e L, 

A e L 6 , 

A G R, 
A G I, 



(4.99) 
(4.100) 

(4.101) 

(4.102) 
(4.103) 
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and 



e 2ite{x) R{\)\ < Ce~ Atx > z < Ce- 4eZ \ XeL 



2„,2 



(4.104) 



where I is an arbitary positive integer and r = £Aq. 



Finally we extend the Riemann-Hilbert problem (13.1 lj) to the aug- 
mente contour Si of Figure 4. From problem (14 .4p and formulae (l4.5p . 
the Riemann-Hilbert problem across K oriented as Figure 3 is given by 

N? = N W (b_)-l s (b + Us 
N w -» I, A -)■ oo. 



(4.105) 



where 



(b±) x ,t,s = 5¥e- u9 ^b±, 
1 fx 



I + w + = < 



I + w. 



1 

1 h 

1 

1 /a 

1 

1 /4 

1 



A < A 
A < A 

A < A 

A < An 



(4.106) 
(4.107) 



(4.108) 



Orient Ei in Figure 4 as in Figure 6 and write 

i (h.y 

1 

1 (hf 



b + = b°b\ = (I+w° + ){I+wl) = < 



1 (h + R) 1 

1 

1 (h 2 + R) 3 

1^0 1 

(4.109a 

1 {h,f \ ( 1 {h 2 + Rf 

1 J I 1 

1 {h,Y \ f 1 (h 2 + RY 

1 J \ 1 

(4.109b) 

where the function {Ky denotes the function is defined in the case j, 
J = 1,2,3,4. 



b_ = b°_b a _ = (I+wZ)(I+w°) = < 



, A < A 

, A > A 
9a) 

, A < A 
, A > A 
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Figure 6. The orient contour of Ex- 

Setting 



i 



' NW\\) =iV«(A), A G 5 U 6, 

=NM(\)(VL)-y s , A G 2 U 4, (4.110) 

=Ar«(A)(^U 5 , AG1U3. 
we find that a simple computation shows that ( 14. 105[) is equivalent to 
the factorization problem 

N?\\) = N^'\\)(J N(1/) (\)) xAS , A 6 Ei, 
iV( 1# )(A)->I, A^oo. 

where 

= (6")^, AgL, (4.112) 

Indeed to show that {b^.)~ tS , for example, converges to I as A — > oo in 
1 we observe that, for fixed ar, t, by formula ( 14. 27ft and the bounded of 
the function 5(A), we have 



\ S * e -*»Wh2(\)\ <c\f3(\)\e- tRaeW 



t 

-oo 



<c^^t \WS)(s)\ds< C 



(4.113) 



|A-i|«+ 2 i-oc lv ' " 71 " |A-i| 
and 



|«V""WJ«A)| < C|S y! i ^ A ° ) ' 1 £ iX^Jji. (4-114) 
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which is converges to as A — > oo, and so on. 
Set 

(^'^ = ±((6^)^-1), (4.115) 

{w 1 ') x , t ,6 = (w+)x,t,8 + {w^) x , t>5 , (4.116) 

Observe from Proposition 14. 71 that . for fixed x,t, we then have 

(u4'W> K'W, e ^(SO n L°°(Si). (4.117) 

4.2. The second step. In this section we show how to convert the 
Riemann- Hilbert problem (I4.11ip on Si to a Riemann-Hilbert problem 
on a truncated contour with controlled error terms. 
From the above section we have 

q(x,t) =2i lim (XN(x,t, X)) 12 

A— >oo 

=i lim A[cr 3 , N(x, t, A)] (4.118) 

A— >-00 v ' 

=i lim A[a 3 ,iV (1) (:r,t,A)] 

A— »oc 

In particular we can take the limit as A — > oo in 5, where N^(x, t, A) = 
N^'\x,t,X), so 

q(x,t) =i lim A[cr 3 , (x, t, A)]. (4.119) 

A— >oo 

The Riemann-Hilbert problem (14.1 lip can be solved as follows (see, 
for example, [15]). Let 

(C±/)(A) = / t/^t^, AeEJe L 2 (E), (4.120) 

denote the Cauchy operator on Ei oriented as in Figure 6. Thus, for 
example, for A > A we have (C+/)(A) = lim e;o J Sl C -{1-^) |^ etc - As 
is well known, the operators C± are bounded from L 2 (Ei) to L 2 (Ei), 
and C + — C- = 1. Also, by scaling, we know that the bounds on 
C± : L 2 (Ex) — )■ L 2 (Ei) are independent of A . 
Define 

G w v J = + (4-121) 
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for 2x2 matrix- valued functions /. By property (14.1 17p . G v is a 

x,t ,5 

bounded map from L 2 (Ei)+L°°(Ei) into L 2 (Ei). Let /i 1 ' = fi 1 ' (X;x,t) G 
L 2 (Ei) + L°°(Ei) be the solution of the basic inverse equation 

^'=I + C wl/ //. (4.122) 

x,t ,5 

Then 

N^(x,t,X)=I+ [ ^ (C;M) 5^ (0 — , AGC\E 1? (4.123) 
y El Q-X 2m 

is the unique solution of the Riemann-Hilbert problem (I4.11ip . Indeed, 

Hi =I + C ± {ii v w 1 l t>& ) 

=1 + C ± {^{wl) x , t , s ) + C ± {^'{w l X^) 

(4 124) 

=I + C wl , fi 1 ' ± fi 1 ' (wl) xAS 

x,t,5 

=// (b±) x ,t,s 



by equation (14.1221) and formula (14. 112ft . which implies that 

N (i') = NP^XX^Us = N?\j N(1/) ) x ,t,s, 
as desired. Substituting formula (14. 123ft into (14.119ft . we learn that 



q(x,t) =- ( / [a 3 , fi 1 ' (C, x,t)wX(C)]— )i2 

J Si 

= - ( / ((i- C7 w r^cK^CC)]^)* 

7 Sj 27U 



(4.125) 



Let w e : Ei — >• M(2, C) be a sum of three terms 

w e = w a + w b + w c . (4.126) 
we then have the following: 

w a = w^ tS \R is supported on M and is composed of terms of type 

w b is supported on L U L and is composed of the contribution to w^ t s 

from terms of type hi- 
w c is supported on L £ U L £ and is composed of the contribution to w^ t S 
from terms of type R. 

(4.127) 
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Set 




Figure 7. The orient contour of cross E^. 

Ei = £i\(MUL £ Ul £ ) (4.128) 

where 

L £ : {A = A + X ue~ l ^;u > e} U {A = A + \ ue l ^~; u > e} 

with the orientation as in Figure 7. Define w' through 

w^ s = w' + w e . (4.129) 

Observe that w' = on S^E^. 

The following estimates are immediate from the bounded of <5(A) 
and Proposition 14.71 in which the decay rate / can be chosen to be 
arbitrarily large. Of course L 2 estimates follow immediately from L 1 
and L°° estimates. However, throughout this article and in particular 
in the lemma that follows, we write out the L 2 estimates explicitly for 
the reader's convenience. 

Lemma 4.8. For A < M, 

H^IU^^nL^Mjai^pR) < Ct~\ (4.130) 

H u ' 6 ||L 00 (LUL)nL 2 (LUL)nL 1 (LUL) < Ct , (4.131) 

lk c ||L-(L e uL e ) < Ce~^\ 

lk c |L 2 (L £ uL £ )<CA|e-^, (4.132) 
\\™ c \\l\u\jU) < C\ e-^ T , 
where j £ = min (4e 2 , ||^). 
Also, 

lki,5llL 2 (Ei) ; IKIUaCEj): Ik'Hi 2 ^) < C. (4.133) 
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Proposition 4.9. In the case, Ao < M and r — > oo, (1 — CV,/) -1 : 
L 2 (Ei) — > L 2 (Si) exists and is uniformly bounded: 

11(1 - C^Witfr) <C. (4.134) 

Corollary 4.10. In the case, Ao < M and r — >■ oo, (1 — C w i> ) _1 : 
L 2 (Ei) — > L 2 (E X ) exists and is uniformly bounded: 

||(1 -C^ r'W^^C. (4.135) 

x,t,o 

A simple computation shows that 

+ / Ei ((l-C 1 ,)- 1 (C 1l ,,I)X 

+/ Ei ((i-c^)- 1 c^(i-c < j- 1 ) 
= f^ai-c^y^w' + i + ii + iii + iv. 

(4.136) 

In the case Ao < M, from Lemma 14.81 it follows that 

\I\ < \\w a \\ L i m + \\w b \\ LHLuL) + |M| L i (LeU £ £ ) 

< Ct~ l + Ct~ l + CAor-' 

< CA r"', as t l > C(t 1 /X ) for A < M, (4.137) 

^ihi^oKJU 2 ^), by mm- 

As above, we have 

IK|| L 2 (Sl) < ct~ l + cT l + cv^r-' < c^/\~ r- 1 , (4.138) 

and 

lk^,5lk 2 (Ei) < IKIU 2 (£ 1 )+||'«/||l2(e 1 ) < c>/\,T- l +c(t\ )-^ < c(t\ )-^(i+A- 1 ). 

(4.139) 

Thus, 

\II\<cX t- {1+ ^(1 + t^- 1 ). (4.140) 
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Continuing, we obtain 

\m\ < IKi-c^^ii^^iic^in^^i^i^^) 

and 

\iv\ < m - crjn^ui - cj )IU 2 (E 1 )ll^ e || L2(Sl )||c . i|| L2(Sl) || 



x ,t ,5 



< c||w e || jL oo (Sl) ||wl' t;(5 || jL 2 (Sl) 



(4.142) 

From Lemma |4"75} for A < M, however, we have ||w e ||i / oo( Sl ) < cr~ l 
and 

lKt,J£a (El) ^ ^IKIli^) + Ikli^sx) + IKIIi^o + lk'll^ (El) ) 

< c(t~ 21 + t~ 21 + AoT-i- + (tAo)-s) < c(*Ao) - i 

(4.143) 

which implies that 

< cr^(tA )-5 (4.144) 

To summarize, for Ao < M, 

\I + II + III + IV\ <c\ t- 1 (4.145) 

as r — > oo. 

We have proven the following result: 

Lemma 4.11. 

q(x,t) = (f ia 3 ,((l-C w ^)(C)wXC)]^hi + O(\ T- 1 ). (4.146) 

We now show that, in general, one may always choose to add or to 
delete a portion of a contour, on which the jump is I, without alter- 
ing the Riemann-Hilbert problem in the operator sense (see identities 
(I4T561) - ( r4~T58l) below). Suppose that £ (1) and £ (2) are two oriented 
skeletons in C with 

card(S (1) n £ (2) ) < oo ; (4.147) 
let u = u(X) = u + (X) + w_(A) be a 2 x 2 matrix-valued function on 
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with entries in L 2 (£< 12 )) n L°°(£ (12) ) and suppose that 

u = on £ (2) . (4.149) 

Let 

# s( i) denote the restriction map L 2 (S (12) ) L 2 (£ (1) ), (4.150) 

I S (i)^ s( i2) denote the embedding L 2 (S (1) ) -> L 2 (£ (12) ), (4.151) 

C* 2 : L 2 (S (12) ) ->■ L 2 (S (12) ) denote the operator in fl4TT2TD with u^w 

(4.152) 

C\ : L 2 (S (1) ) L 2 (S (1) ) denote the operator in fl4TT2TD with u t <R> w*^, 

(4.153) 

: L 2 (£«) -> L 2 (S< 12 )) denote the restriction of C* 2 to L 2 (£«). 

(4.154) 

And, finally, let 



i' 

x,t,8i 



(4.155) 



I S (i) and I S (i2) denote the identity operators on 
L 2 (£«) and L 2 (£( 12 )), respectively. 

We then have the next lemma: 
Lemma 4.12. 

Cf = C E U Cf , (4.156) 

(I E (i) - C^y 1 = i2 s( i) (I E (i2) - C^ 2 )~ 1 Ii;(i)^s( 12 ) ) (4.157) 

(I s(12) - ClY 1 = I s(12) + (I s(1) - Ci)- 1 ^,,), (4.158) 

m t/ie sense £/ia£ if the right-hand side of ( 4-15T\) ,resp. exists, 



then the left-hand side exists and identity \J^.151\j ,res'p. \J^.159^ , holds 
true. 

Proof. The proof of identity (I4.156P is trivial. If g E L 2 (S^ 12 ^) and 
A G then 

((I EW - C l u )R^ )g )(X) = g(X) - (ClR^gKX) 

= g(X) - (Cl*g)(X) (4.159) 
= {{l^-C^)g){\). 
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Hence, for / 6 L 2 (S^^), we have 

((I E( i) - C^)i2 E (i)(I E (i 2 ) - C^ 2 ) _1 I S (i)^ s( i2)/)(A) = ((I S (i2) - C^ 2 ) 

((I E(12 )-^ 2 )- 1 I s(1) ^ s(12) /))(A) 

= (I E (i)_ E ( M )/)(A) = /(A). 
(4.160) 

Conversely, if / e L 2 (£ {1) ) and A e £W, then 

(I S (l)^ S (12)(I S (l) - C u )f)(X) = ((I S (1)^ E (12)/) - C^ 2 (I S (1)^ S (12)/))(A) 

and hence, 

(Is(i)-».S(l 2 ) Ot(l) _ Cu)/)(A) = (Is( 12 ) - C« 2 )(Is( 1 )^S( 12 )/) + #2 

where g 2 = on However (I S (i2) — C]?)g 2 = I E (i2)<72 = g 2 , and so 

(I E (i2) — C^ 2 y l g 2 = g 2 . It follows that 

J? E (i) (I S (i2) - C^ 2 ) _1 I E (i)^ S (i2)/ = (I s( i2) - C'^ 2 )~ 1 (((I S (i2) - C,l 2 )Is(i)^i;( 12 )/) + #2) 

= i? S (i)I S (i)^ E (i2)/ + R Ti d)g 2 = f. 

This proves identity (14.1571) . 

On the other hand, using (I4.156p . we get 

(I s(1 2) - 0(I 2(1 2) + C^(I EW - C^R^ W ) 

= (I E (i2) - C* 2 ) + C u (I E (i) - C^)(I S (i) - C^) _1 i? s( i) 

= I E (12) - C^ 2 + Cf i? S (l) = I S (12) 

Conversely, i? E (i) (I S (i2) — C^ 2 ) = (I S (i) — C^)i2 E (i), and so we have 

(I 2( i2) + C u (I S (i) - C^) _1 -R E (i))(I S (i2) - C* 2 ) 

= (I E (i2) - C* 2 ) + C u (I 2 (i) - C^)(I S (i) - C^) _1 i? s( i) 

= I 2 (12) - C^ 2 + Cf i? E (l) = I S (12) 

This proves identity (14.1581) and the lemma. □ 

We apply Lemma 14.121 to the case £« = S' 1; S( 12 ) =Ei and u = w . 
We learn in particular that 

the boundedness of ||(I — C w i)~ l ||l 2 (£i) is equivalent to ^ 
the boundedness of ||(I — C' w ,)~ l || L 2( 2 ' ) 
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(Here I' = Ijv , C' w , = C w i .) Also, as in the proof of Lemma 14.121 from 
identity (14. 157ft we have 

(I'-CV)- 1 =Rv i (I-C w ,y\I + g 2 ), g 2 = 0onEi 

= R^(I - C^)- 1 ! + ife^ (4.162) 
= R^(I-C wl )-H. 

Inserting this identity in formula (I4.146P yields the following proposi- 
tion: 



Proposition 4.13. 

?(*,*) = (/ [t73,((l / -C;,)- 1 I)(CV(C)]^)2i + O(A r-'). (4.163) 
Set 



E / 27ri ' 



L' = L\L e . (4.164) 

Then 

S; =L'U 11 ■ (4.165) 
On S; set // = (V - C^) -1 !. As in formula (143231) . it follows that 



solves the Riemann-Hilbert problem 

ivf ) (A) = iv: (1 ' ) (A)(j; (1 , ) )^(A), ags; 

Af' (1 ') (A) ->■ I A^oo. 



(4.166) 



(4.167) 



where 



u/ = u4 + u;l, (4.168) 

&' ± = ]I + u4, (4.169) 

(&±W = ^ 3 e-^6' ± , (4.170) 

(^)).A* = (0^+W (4-171) 
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From formulae (14. 109[) and ( 14. 112ft we can get that 
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b' 



1 (RY 

1 

1 (Rf 
1 



A < A 
A > A 



V 



1 
1 



b' 



1 
1 



b' 



1 (i?) 2 

1 

1 (i?) 4 

1 



A < A 

A > A 
(4.172) 



4.3. The next step. Extend S' x to the contour 
S; = {A = Ao + Aowe^w G M}U{A = Xo+Xoue^, u e R}. (4.173) 
Define the scaling operator 



/(A) h> Nf(X) = f{—7= + Ao) 

vol 



Then we can get 
q(x,t) = 



(8t)- i "e mx ° +2l ~<(n° 1 ) 12 + o(r 



(4.174) 



(4.175) 



where 7 = ± j*l log |A - A |dlog(l - A|p(A)| 2 ). 

We notice that in the neighborhood of the stationary point A = Ao, 
the function 5(A) appearing in the formula (14. 2p can be represented 

as(H) 



5(X) ± = (A - Ao)^^^ 1 ^^'^^ 10 ^ 1 ^'^')! 2 ) (4.176) 



where 



v = ~— log(l-A |p(A )| 2 ) (4.177) 

and (A — Ao)± denotes the boundary values of the corresponding mul- 
tivalued function defined on the A— plane with the cut along (— oo, A]. 
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Then a straightforward computation shows that as t — > oo 

{N5 a3 e- itea3 [f]){\) ^<f) &3 \ uia3 e- l ^ aa [f}{\ ), (4.178) 

where 

(f)= (8t)~T e 2 ^o e ^/-Sc lo s( A o-C)rfiog(i-C|p(C)l 2 ) (4.179) 

and [/](Ao) is defined by [f] = R + hi where R and h\ are defined in 
(14.11) in subsection 14.11 

It follows from the exponential decay of e ~ CT3 [/](Ao), that the as- 
ymptotic formula in (I4.178P has an L l PI L 2 PI L°°(Si — Ao) error of 
order Since is indenpendent of A, iV° is the solution of the 

Riemann-Hilbert problem on Si — Ao, 

N° ->• I, as A ->■ oo. 

if and only if <p aA N Q is the solution of the Riemann-Hilbert problem for 
the jump matrix given by the right-hand side of (I4.178p . Deforming 
the Riemann-Hilbert problem (14. 18QH on Si — Ao to the real axis we 
obtain the Riemann-Hilbert problem 

N° -> I as A -> oo. 

which can be solved in closed form. This problem was first considered 
by Its, and the following calculations can be found in [7] [8] [6]. 
Setting 

V>(A) = Ar (A)A^ ff3 e-^ ff3 , (4.182) 
we can represent the Riemann-Hilbert problem (I4.18ip as 

^ + (A)=^_(A)[/](A ), AgK, 

,2 (4. led J 

V»(A) ->■ A^e"— CT3 , as A ^ oo. 

By differentiation we have 

( ^ + ^ W)+ = ( ^ + ^ W) - [/](Ao) ' AGM ' (4 - 184) 
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Now as det [/](Ao) = 1, it follows by a Liouville argument that det if) = 

■it _l ±a 

d\ ~ 2 l 



1. Hence, if) 1 exists and is bounded. But (4f + \iXa 3 if))if) 1 has no 



jump across K. and must be entire. Also 

= OiX-^ + ^N?] 



It follows by Liouville 's argument that 

dip 1 
dX + 2 



(4.185) 



<IU + \iXazif) = flip, (4.186) 



where 



In particular, 



(jV°) 12 = _z/3 12 . (4.188) 
Consider first ImA > 0. From equation (14.1861) we obtain 

d 2 ih + A 2 ? 

^ = (-^-£ + /WWfi- (4-189) 

Setting 

^(A) = s(e"*A) (4.190) 

results in the reduction of equation (14.1901) to the Weber's equation 

|l + - J + °M0 = 0, (4-191) 

where a = ifiufoi- 
Hence, 

Vi(A) = ClJ D a (e-^A) + c 2J D a (-e"^A) (4.192) 

where D a {-) denotes the denotes the standard (entire) parabolic-cylinder 
function. 
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From [IT] , we know that as ( — > oo, 
A,(C) = C a e~^ 2 (l + 0(r 2 )), I argd < f , 

= Ce"K 2 (l + 0(r 2 )) - i^e^C-^e^'Cl + 0(C 2 )), f < argC < f , 

= C a e"^ 2 (l + 0(r 2 )) - j^e-^C-^^^^l + 0(r 2 )), -f < argC < -f , 

(4.193) 

where T is the Gamma function. 

Setting A = e^V, with a > 0, comparing the right-hand side and 
the left-hand side of (14. 192ft . we conclude that 

37VV 

c 2 = 0, a = iv and C\ — e 4 , (4.194) 

so that 

^+(A) = e-^Daie'^X), ImA > 0. (4.195) 
From equation ( 14.1861) we learn that 

= ^ e ~^(J^(e-^A)) + |zAD a ( e -^A)),ImA > 0. 

(4.196) 

In a similar way we find that 

^(A) = eTD_ a ( e -TA), ImA>0, (4.197) 

and 

^i + 2 (A) = -^e^(A (jD _ a(e -f A)) _ W a (e-TA)),ImA > 0. 

(4.198) 

Repeating these calculations for ImA < yields 

^h(A) =eTD a ( e fA), (4.199) 
«A) = ^ e¥ (^(^( ef A )) + iiAU.(e*A)), (4.200) 
V&(A) = e-^L»_ a (e^A), (4.201) 
lfo(A) = ±e-^^(D. a (e^X)) ~ i<AD_ a (e^A)X4.202) 
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Substituting these formulae in (14. 183[) gives 



W_)-V + = /(A„)= ('-^f ' P \ X0) ) (4.203) 

V A p(A ) 1 J 

then we can obtain 

Aop(Ao) = tpuipti ~ fax&ti 

= (e^D a (efX))(^e-^(^(D a (e-^X)) + \r\D a {e~^ X))) 

-^(i( D ^X)) + ±i\D a (e^\))e- 3 -¥D a (e-^\) 
= ^W(D a (ef\)),D a (e- 3 -?\)) 

~~ /9i 2 r(-o) • 

(4.204) 

where W(f, g) = fg' — f'g is the Wronskian of / and g. 
Thus 

ymer^e 2 

/5i2 = , 7T , v( S 4.205 
A p(A )r(-a) 

Finally, careful bookkeeping of the error terms and the Lemma 8.1 
in [19] yields the following result. 



Theorem 4.14. Letq(x,t) be the solution of the Cauchy problem M.l\) . 
Then as t — >• 00 

q(x,t) = q as (x,t) + 0( 1 -^) (4.206) 



where 

qas = ^ t a(Xo)e^-^ x ^\ 

|«(Ao)| 2 = ^ = -^log(l-A |p(A )| 2 ), 

arga(Ao) = -3z/log2 - f + argT(^) - argr(A ) + J log |A - A |dlog(l - A|p(A)| 2 ), 



\ — x 



(4.207) 
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